dt=0.01;
t=0:dt:3;
f1= heaviside(t) - heaviside(t-2);
f2 = heaviside(t) - heaviside(t-3) +  heaviside(t-1) - heaviside(t-2);
f=conv(f1,f2)*dt;
n = length(f);
tt = (0:n-1) * dt - 2;
subplot(221);
plot(t,f1);title("f1=heaviside(t)-heaviside(t-2)");xlabel('t');
subplot(222);
plot(t,f2);title("f2=heaviside(t)-heaviside(t-3)+ heaviside(t-1)-heaviside(t-2)");xlabel('t');
subplot(223);
plot(tt,f);title("f = conv(y1,y2)*dt");xlabel('t');